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Surfaces of surfactant-templated polyacrylamide hydrogels were imaged by atomic force mi- 
croscopy (AFM), and the surface morphology was studied by numerical scaling analysis. The tem- 
plated gels were formed by polymerizing acrylamide plus a cross-linker in the presence of surfactants, 
which were then removed by soaking in distilled water. Gels formed in the presence of over 20% 
surfactant (by weight) formed clear, but became opaque upon removal of the surfactants. Untem- 
plated gels formed and remained clear. The surface morphology of the gels was studied by several 
one- and two-dimensional numerical scaling methods. The surfaces were found to be self-amne on 
short length scales, with a roughness (Hurst) exponent in the range 0.85 to 1, crossing over to a 
constant root-mean-square surface width on long scales. Both the crossover length between these 
two regimes and the saturation value of the surface width increased significantly with increasing 
surfactant concentration, coincident with the increase in opacity. We propose that the changes in 
the surface morphology are due to a percolation transition in the system of voids formed upon 
removal of the surfactants from the bulk. 

PACS numbers: 82.70.Gg, 68.37.Ps, 61.43.Hv, 05.40.-a 

I. INTRODUCTION 

Polyacrylamide (PAAm) gels are used extensively for separations of biological macro-molecules. Because of their 
intrinsically broad distribution of pore sizes they are very useful for a wide range of separations for protein isolation 
and sequencing of single-stranded DNA jlj . In particular, protein separations achieved by electrophoresis on gels tem- 
plated with surfactant micelles have recently been reported |^] , and the possibility of improvements of the separation 
properties, when pores of a particular size are templated within the gels, have been proposed Extensive work has 
also been done on the possibility of using such templated gel matrices in gel permeation chromatography (GPC) Q. 

Since its introduction in 1986 ||, atomic force microscopy (AFM) has been widely used for recording images of many 
surfaces with resolution at small length scales. Recent developments in AFM technology have enabled its use, not 
only in high-resolution profiling of surface morphology and nanostructure, but also in determining local mechanical 
properties and in surface compositional mapping of heterogeneous samples. These applications have been extensively 
discussed by Magonov et al. || Despite such developments, there have been only a few reports in the literature 
on the direct observation of polymer gel surfaces by AFM M, 0, |l(J. Here, an AFM study of the surface morphologies 
of templated PAAm hydrogels and their scaling properties are presented. 

We use micelles of tetradecyltrimethyl ammonium bromide (TTAB) surfactant to template pores in PAAm. Aque- 
ous solutions of surfactant, acrylamide, and N',N'-methylene bisacrylamide (cross-linker) are mixed at various con- 
centrations that allow the surfactant molecules to self-assemble into micelles. The solutions are then polymerized, 
after which the surfactants are removed by soaking in distilled water to produce templated hydrogels. A schematic 
representation of the templating process is shown in Fig. |l|. In previous applications of TTAB templated PAAm 
gels for GPC Q, mean pore sizes, exclusion limits, and mechanical integrity were determined. It was shown that 
templated gels yield improved resolution for the separation of globular proteins from 10,000 to 150,000 Da molecular 
weight. Electrophoretic separations achieved on templated PAAm were also shown to be consistent with the pore-size 
distributions observed by GPC analysis. Templated PAAm gel matrices have the potential to improve separation 
efficiency over a wide range of globular protein sizes. 

In this paper we present a quantitative study of the structural changes due to the templating process, using scaling 
analysis of AFM images of the surfaces of templated and untcmplated gels. A wide variety of image- analysis tools 
are available with the Digital Instruments AFM software |ll|. However, these are not well suited for scaling analysis 
of soft surfaces. Hence, we investigated the surface structure over a wide range of TTAB concentrations using our 
own, custom-designed surface-analysis software. This scaling analysis appears to be a robust tool to analyze and 
characterize surfaces imaged by AFM. 

The rest of this paper is organized as follows. In Sec. II we describe the experimental procedure, apparatu s, and 



results. In Sec. Ill we disc uss th e scaling-analysis methods. This section is subdivided into two-dimensional (Sec. [II A ) 



and one-dimensional (Sec. IIIB| ) methods. In Sec. [V we discuss the results of the scaling analysis of the AFM images 
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and in Sec. [V| we present our conclusions. 

II. EXPERIMENTAL 

PAAm gels were formed at room temperature by free-radical polymerization with chemical initiators from a mixture 
of monomer, cross-linker, buffer, and surfactant. The gels were cast between two glass slides separated by 1.5 mm 
spacers. The glass slides were ordinary precleaned microscope slides (Fisher Scientific). They were cleaned with 
detergent and rinsed in distilled water and air dried before use. AFM images (not shown here) of these slides revealed 
an root-mean-square (rms) roughness of less than a nanometer over several square micrometers. Acrylamide (Elec- 
trophoresis Grade, Fisher Scientific) solutions were polymerized in the presence of a TBE buffer (45 mM tris-borate, 
1 mM EDTA) at pH 8.3. The chemical initiators were ammonium persulfate (Ammonium Peroxydisulfate, Certified 
A.C.S., Fisher Scientific) and N, N, N', N'-tetramethylethylenediamine (TEMED, Fisher Scientific). The concentra- 
tion of acrylamide was kept constant at 40% by weight, and that of the cross-linker, N',N'-methylene bisacrylamide 
(Bis-acrylamide, Electrophoresis Grade, Fisher Scientific) was held at 7% of the acrylamide concentration in all sam- 
ples. Various amounts of the surfactant, TTAB (Sigma- Aldrich), were added to form the templated gels. The amount 
of surfactant added was so that the final surfactant concentration ranged from 0% to 40% by weight in steps of 10%. 
The presence of surfactant micelles did not prevent the formation of PAAm gels. The rates of polymerization were, 
however, altered by the presence of surfactants. Higher concentration of surfactants in the pre-gel solution reduced the 
rate of acrylamide polymerization. All pre-gel solutions were clear, homogeneous solutions with viscosities dependent 
on the amount of surfactant. The resulting templated gels were all optically clear before surfactant removal. 

After polymerization, the gels were peeled off the glass plates and soaked in distilled water to remove the surfactants 
by diffusion. More than 98% of the surfactants diffused out of the gel upon soaking for about two days. Several 
methods, including Raman spectroscopy have been used to quantify the removal of the surfactant templates by free 
diffusion into water Q. During the surfactant-removal process, gels with surfactant concentrations of more than 20% 
transformed into white homogeneous materials. This suggests a structural change in the templated gels, involving 
a length scale comparable to the wavelength of visible light. This further indicates that the structural change is 
associated with kinetic restructuring of the gel during and after the removal of micelles. 

The gel surfaces were imaged using AFM (Digital Instruments, D3000) with a silicon nitride probe (Digital In- 
struments) of spring constant 0.32 N/m. Scans were performed in contact mode under water (HPLC grade). The 
set-point voltage was kept at the lowest possible value to minimize the risk of damage to the gel. The inherent soft 
nature of the gels under investigation can cause difficulties in imaging and affect the reproducibility of the images. In 
order to overcome this difficulty, we devised a sample holder which helped provide the required stability and hence 
improved the reproducibility of the images. This sample holder consists of a stainless-steel block with a well in the 
center, into which the sample is placed. The sample is then covered using a thin steel cover slip with a hole in the 
center for imaging. The cover slip is held in place with mounting screws (see Fig. |^). 

After the surfactants were removed by soaking, three different locations on each gel surface were selected for imaging, 
and a set of centered, zoomed images were obtained for each location. The image sets were created by scanning the 
surface from small scan sizes up to larger sizes, in sequence. The images were all L x L squares with sides L=1.25, 
2.5, 5.0, 10.0, and 20.0 /im, respectively. This process yielded a set of successively zoomed-out images of the same 
location. The lateral drift of the sample was very small. Figure || shows one such zoomed set, and the excellent 
reproducibility of the images justifies the use of the cover slip holder. To further test the reproducibility, drift, and 
surface distortion, two images were obtained for each scan size, first in the forward scanning direction, and next in 
the reverse scanning direction, before changing to the next scan size (see Fig. ^). Some distortion along the direction 
of the tip motion is noticeable. All images have 512 x 512 pixel resolution, regardless of scan size, and all imaging 
was performed at room temperature. 

This telescoping scanning procedure produced a set of ten images for each location and a total of thirty images for 
each gel. With this procedure, five different gels were examined with surfactant concentrations of 0, 10, 20, 30, and 
40%, yielding a total set of 150 images (75 forward scans and 75 reverse scans). Many of the images were too noisy 
for analysis, and all 150 images were examined for quality in a blind procedure in which two of the authors, who had 
not themselves acquired the data, were separately shown all 150 images. They indicated whether each image should 
be retained for analysis. If an image was indicated as too noisy by either individual, that image was excluded from 
analysis. This procedure yielded a total of 106 images for analysis, as detailed in Table || 

The gel surfaces were also imaged using tapping mode under water. A single set of ten images as described above 
was obtained for gels with surfactant concentrations of 0% and 40%. The surface morphologies were statistically 
indistinguishable from those obtained using contact mode. Since tapping mode is rather time consuming and did not 
provide any apparent advantage over contact mode, all the data presented here were collected in contact mode. 

Figure shows examples of AFM images (L = 10.0 fim) of gel surfaces with each of the five different template 
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concentrations. The gels exhibit marked statistical differences, most notably that the height range of the surface 
depends on both the surfactant concentration and the length scale of the scan. This is clearly seen in the three- 
dimensional visualizations of the gel surface, shown in Fig. ^|. Typical scan lines for gel surfaces at various surfactant 
concentrations are shown in Fig. |7J. On the 10 /im length scale, the untemplated gels have height ranges on the order 
of a few nanometers, while the 40% templated gels have height ranges on the order of a few hundred nanometers. 
AFM images and Raman spectroscopy of glass plates before and after polymerization and gel removal indicated that 
there was no polymer left behind on the glass. This suggests that the gel surfaces imaged by AFM are relaxed, free 
surfaces, rather than fracture surfaces. 

The surface analysis routines available with the D3000 only quantify changes in the overall surface width (often 
called "the surface roughness"). In contrast, our statistical analysis shows self-affine scaling behavior at small length 
scales, which crosses over to a length-scale independent behavior at large length scales, as discussed in Sec |fv|. As 
we are interested in the behavior of the surface on a wide range of length scales, no image filtering or processing was 
performed, with the exception of flattening (described in Sec. [II A). It is our opinion that scaling analysis of filtered 
images could be highly suspect since the essence of most image processing techniques is to alter the data differently 
on different length scales. All analysis was performed with our own code, written in C, as discussed in Sec. III. 



III. ANALYSIS METHODS 



We here present a well-defined procedure to characterize surfaces, based on the statistical analysis methods described 
m Refs. H One common measure of surface morphology is the rms height of the surface, which is sometimes 
referred to as the "rms roughness" (referred to here as the rms surface width) , 



w 



\ 





Ei=i h i 


M 2 


M 2 



(1) 



where hi is the ith value of the height in an image consisting of M x M pixels. (For all of our images, M = 512.) 
However, such a drastic simplification of the surface structure to only one single number ignores nearly all of the 
morphological information contained in the image. Reporting the surface width without specifying further details can 
be meaningless. The value of the width is highly dependent, not only on the image processing methods used |fL4| , but 
also on the length scale of observation ]l5| . We therefore present a comprehensive scaling analysis of the gel surfaces 
which goes beyond the simple single w description. 

The surface width defined by Eq. ([[]) for a L x L (fim 2 ) image of M x M pixels varies with L. For a variety of 
different surfaces the width is found to exhibit different length dependences with two limits: 
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where the exponent a is the Hurst exponent (in this study equal to the roughness exponent Jl6| ), w sat is a constant 
saturated surface width, and l x is a cross-over length scale between the two behaviors. In general, a, l x , and w sa t 
depend on the details of the surface, such as the surfactant concentration in the present study. 

The power law in Eq. (j^) is a consequence of self-affine scaling behavior, in which the surface is statistically invariant 
under the transformation fl2|| 



bx, y — > by, h — > b a h 



(3) 



where x and y are distances in the two image-scan directions, h is the height, and b is a dimensionless re-scaling factor. 
Self-affine behavior is different from self-similar behavior, which involves invariance under the isotropic scale-change 
transformation ]12] | x — > bx, y — > by, and h — > bh. This distinction between self-affine and self-similar structures is 
often overlooked, and the term "fractal" is used somewhat freely for both in the literature. 

If the surface were ideally self-affine, then the power-law behavior in Eq. (^) would continue on all length scales, 
and there would be neither a saturation width, w sa ±, nor any identifiable length scale such as l x or a typical feature 
size. For a surface whose heights are described by an ideal random walk, a = 1/2 and l x —> oo. For most surfaces, 
measured values of a are between 1/2 and 1 JIt], [ll| O], even for such different surfaces as brittle rocks (a « 0.78) |2C(| 
and CuCl islands on CaF2 (a « 0.84) ]2~]] ]. As discussed in Sec. IV, our estimates of a for the different gel surfaces 
all fall within the range [0.5, 1.0], and w sa ± and l x vary strongly with the surfactant concentration. 

In the following subsections, we discuss several different methods used to examine the scaling behavior of the gel 
surfaces. Both two-dimensional analysis methods that analyze the image as a whole, and one-dimensional methods 
that analyze the image one scan line at a time were examined. 
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A. Two-Dimensional Methods 



The AFM records the vertical surface height, h(x,y), as a function of the fast (x) and slow (y) scanning directions. 
As already discussed, L dependent surface widths can be obtained from L x L (/mi 2 ) images. After imaging at a 
sufficiently large number of scan sizes (different L's), a, w sat , and l x can be obtained. However, the number of images 
required for such an analysis is large. By noting that each image contains information about length scales from L/M 
(where M is the number of pixels) up to L, one can perform analyses over a wide range of length scales from a single 
AFM image. Several variations of two-dimensional analysis methods were studied. We discuss a spectral method 
based on the Fourier transform of the image Jl3|], a box-counting method from subsections of images, and widths 
calculated from the same subsections of the image |L2] . All of the images were flattened using the AFM software |ll[] 
and not processed any further. The flattening procedure fits a straight line to each scan line by the method of least 
squares, which is then subtracted from the data. The fitting and subtracting procedure was applied, independently, to 
each scan line in the image. After flattening, the average height and slope in the x-direction (the fast scan direction) 
are zero. 

For the spectral methods, we computed the structure factor, S(k) oc | T^{h{x, y))\ 2 , which is the absolute square 
magnitude of the complex Fourier transform of the surface height. Here, k is the wave vector in the x, y plane. For 
large k, the structure factor should scale as S(k) cx fc~( 2 + 2Q ) [pi| ; a generalized Porod's Law J22[ ]. The presence of 
anisotropy in the AFM images produces ridges along the k x and k y directions as shown in Fig. gTThis made circular 
averaging unreliable, and hence the results from spectral methods were suspicious. The presence of the ridges seems 
to be a general feature of AFM imaging, as has been observed for hard materials as well jlj]. Removal of these ridges 
(which have a non-zero width) by a technique such as that suggested in Ref. p4| would require a priori knowledge of 
their functional form. 

For both the box-counting and width methods, the L x L (/im 2 ) image of M x M height values (pixels) was 
recursively subdivided into non-overlapping images of I x I (/mi 2 ) with (M/n) x (M/n) pixels, where I — L/n /im 
and n < M/2 is an integer denoting the level of subdivision. For each value of I, (L/l) 2 different sub-images were 
obtained. Non overlapping sub-images were used to ensure statistical independence. 

A box-counting method was tried following Ref. ||l2f . The minimum number of cubic boxes of side I needed to span 
the height range of a sub-image of I x I /im 2 was calculated. However, the height range was nearly always less than I, 
yielding the total number of boxes to be one, regardless of the length scale. This would seem to imply a = 1. Thus 
the box-counting method is entirely insensitive to the scaling behavior of the gel surfaces. 

The rms width for the j-th sub- image of size i x i, 
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and its average, (wi), over all sub-images of size / were calculated. Typical results are shown in Fig. ^. For a discussion 
of the effects of changing the order of the averaging and the square root, see Ref. p3| . 

Among the several two-dimensional analysis methods used to study the scaling properties of the gel surfaces, the 
calculation of (wi) was found to provide the most reliable results. The other methods have inherent disadvantages as 
discussed in the preceding paragraphs. As we shall explore in the next subsection, the scanning process introduces 
noise between the scan lines, and for a more accurate analysis, one-dimensional methods should be used. 



B. One-Dimensional Methods 



If there is noise in the measurement process we would expect to measure h measme d(x,y) = h(x,y) + h no { sc , where 
the noise adds an effective offset. If we further assume that h no i so is a random walk in time, noise correlations in 
time would appear as correlations in space. Furthermore, if the noise varies sufficiently fast compared to the scanning 
velocity, h no i se at the beginning of each scan line will be completely uncorrelated. Thus for a perfectly flat surface 
with h(x, y) = 0, each scan line would appear to be a random walk with a — 1/2, with little or no correlation between 
scan lines. Motivated by these considerations and the extremely small height ranges in the gels with 0% and 10% 
surfactant, we constructed a one-dimensional analysis procedure which analyzes the images one scan line at a time, 
where a scan line is a line in the fast (x) scan direction. 

One-dimensional surface widths (wi) were calculated in a manner equivalent to that defined in Eq. (0) with j now 
representing a sub- interval in the fast scan direction instead of a sub- image and (M/n) 2 replaced by (M/n). The toy, 
for all j and all scan lines, were then averaged to get the (wi) for the image. The surface width calculations at the 
smallest image sub-interval depends strongly on the definition of wi (see Eq. (Eh). When applied to a straight line 
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(a = 1), the standard one-dimensional estimate of (wi) used here overestimates a by as much as 20% for sub-intervals 
that contain less than four data points. The results of the one-dimensional surface width analysis are presented in 
Table || and in Sec. |Fv| 



While surface- width calculations contain information about all length scales < L, calculation of correlation functions 
yield information at a specific value of relative separation. We define the height-height correlation function of a scan 
line at fixed y as 

C(r, y) = (h(x, y)h(x + r, y)) x - (h(x, y)) 2 x , (5) 

where r is a relative displacement in the x direction, and (-^ denotes an average over all values of x in the scan 
line. Note that after flattening, the second term is zero, and that C(0,y) is the rms width of the scan line, w(y). 
Although the flattening was necessary to compensate for non- horizontal mounting of the gel and/or gradients in the 
gel thickness, it introduces correlations at the length scale L. Consequently, correlation functions for flattened images 
are meaningful only for lengths, I < L/2. For scaling analysis, it is more convenient to deal with an alternative form of 
the correlation function, which is sometimes referred to as the variance of increments pi[ or the increment correlation 
function, 

H(r,y) = (\h(x + r,y)-h(x,y)\") x . (6) 

For the special case of q = 2, 

H(r, y) = ([h(x + r,y)- h(x, y)] 2 ) x = 2[w(y) 2 - C(r, y)] . (7) 

The correlation functions in the slow (y) scan direction, C(x,r) and H(x,r), can be analogously defined. 

To improve the statistical accuracy of the increment correlation function measurements, multiple scan lines must 
be used. Thus, we define the first and second moments as 
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fl-fiurt(r) = (H(r, y)) = £* =1 [£^i H(r, y = iL/M)\ /pM , 
(H(r,y) 2 ) = EU [££i H(r,y = zL/M) 2 ] /pM , 

where Y^j=i runs over the p different scan lines. The standard error of (H(r,y)) is then <7f as t = <^(f,y) = 
(H(r,y) j-(H(r,y)) _ Again, the moments and error in the slow scan direction can be analogously defined. 
The increment correlation function for q — 2 is expected to scale as 

Mr)«(2. 15! X ■ (9) 
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In Fig. [l(] we see power-law scaling for small r and an apparent cross-over to saturated values at large r. When the 
noise-induced offsets between scan lines are small, H s \ ow (r) should have the same scaling as -f/fastM- However, if the 
signal-to-noise ratio is small, we expect H s \ ow (r) to measure only the noise with a = 1/2. Figure [ll] shows plots of 
H s i OVJ (r) and H{ ast (r) for 0% and 40% templated gel surfaces. For the untemplated gel surfaces the fluctuations due 
to noise in the slow scan direction are comparable to the signal. Thus H s \ ov/ (r) does not follow Hf ast (r). For the 
40% templated gel surfaces, the random offsets between scan lines are relatively smaller, but not negligible. Hence 
-HsiowM is much closer to H faBt (r). 

By evaluating the increment correlation function and plotting H{ as t(r) on a log- log scale, the roughness exponent 
a, the limiting surface width u> sa t, and the lateral cross-over length l x can be determined from the slope in the 
small-r linear regime (see Fig. |l0|) , the limiting value of Hf ast (r) for large r, and the inflection point of the logarithmic 
derivative of Hf ast (r) (see Fig. |12|), respectively. 



IV. RESULTS AND DISCUSSION 



As shown in the previous sections, among the two-dimensional analysis methods, only the Z-dependent surface 
width (wi) gives meaningful results for this system. The /-dependence of this quantity was presented in Fig. || and 
is essentially similar to that of the one-dimensional surface-width (not shown) and the one-dimensional increment 
correlation function shown in Fig. |l0|. All three quantities show a self-affine scaling region on short length scales, and 
a saturated plateau on large length scales. These two scaling regimes are separated by the crossover length l x . Both 
the crossover length and the saturation value increase with increasing surfactant concentration. 
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The values of the roughness exponent a extracted from the self-affine regime for the surface width measurements 
(both two-dimensional and one-dimensional) , along with those calculated from the one-dimensional increment corre- 
lation function with q = 2, are listed in Table ||. Using q ^ 2 in Eq. (0) yielded the same values of a as obtained with 
q = 2. This indicates that the gel surfaces are perfectly self-affine on these short length scales, rather than multi-affine 
]l6| . |23| , p5| . A plot of the estimated values of a vs the surfactant concentration is shown in Fig. [l3|. The estimates 
obtained from the two-dimensional surface width measurements and from the one-dimensional increment correlation 
function agree well for 20%, 30%, and 40% templated surfaces. We believe that the lack of agreement for 0% and 10% 
templated surfaces is associated with the low signal-to-noise ratio in the images for these gels. The vertical resolution 
limit of the AFM (> 0.1 nm |ll[]) contributes to the uncertainty since the height variations are only about 1 nm for 
the lowest surfactant concentrations. 

The values for a, determined from the one-dimensional surface- wid th me asurements, are consistently higher than 
those obtained from the other two measurements. As discussed in Sec. IIIB, the roughness exponents obtained using 
the one-dimensional surface width for intervals with only a few data points are consistently too high, and this may 
be at least part of the explanation for the discrepancy. The values of a for 20%, 30%, and 40% templated surfaces 
obtained from two-dimensional surface-width and one-dimensional correlation-function measurements are clustered 
between 0.80 and 0.89. However, measurements of surface-roughness exponents from relatively small data sets, such as 
the present ones, tend to have systematic errors, as well as random ones, and the calibration curves for synthetic data 



sampled over 512 points, given in Ref. |26 , indicate that the values obtained here may be underestimated by up to 
10%. Thus, based on our results it is not unreasonable to conclude that all of these surfaces have the same roughness 
exponent in the range 0.85 to 1.0, irrespective of the surfactant concentration. A roughness exponent in this range is 
consistent with previous results for N-isopropylacrylamide gel surfaces by Suzuki et al. [^). However, as discussed in 
the next paragraph, the magnitude of the height variations depends strongly on the surfactant concentration, as do 
the crossover length l x and the saturation surface width w sa t- 

The limiting value of the surface width on large length scales, w sa t, is found to increase by almost two orders of 
magnitude as the surfactant concentration is increased from 0% to 40%. This behavior is evident from Figs. 0, Ufa), 
0(a), and|l|. The measurements of w sa t are listed in Table || and displayed vs surfactant concentration in Fig. 
Similarly, but less spectacularly, l x increases by about a factor of two, as shown in Table |j| and Fig. [l5|. 

As the surfactants are removed, the transparency of some of the gels change. The gels that contained no or 
little surfactant stay clear. On the other hand, gels that had higher concentrations of surfactant (more than 20%) 
become opaque and white as the surfactant diffuses out. Using the theory of light scattering by refractive index 
fluctuations, an expression for the spectral density of light scattered can be obtained. Reasonable functional forms for 
the correlation function can then be used to obtain the turbidity or inverse attenuation length. A detailed derivation 
is provided in the Appendix. The turbidity for a wavelength A = 550 nm in the visible region is shown as a function 
of the correlation length £ in Fig. |l^. The two curves presented correspond to an isotropic exponential and Gaussian 
correlation function, respectively. If we associate the correlation length £ with the saturation value of the rms surface 
width, w sat , we notice that the turbidity increases by about six orders of magnitude as w sa t changes from 1 nm for the 
0% templated gel to 80 nm for the 40% templated gel. These changes in the turbidity are consistent with a change 
from a transparent to an opaque material. Our justification for the association between £ and u> sa t is discussed below. 

The experimental results presented here indicate that while the templated gel surfaces are self-affine on small to 
moderate length scales, with a in the range 0.85 to 1, both the crossover length and the saturation value of the 
surface width increase significantly with the bulk concentration of surfactant. From Figs. |lj and [l^ it is clear that 
most of these increases occur only after the surfactant concentration exceeds about 20%. This indicates that the 
dramatic changes seen in the surface structure are reflections of a sharp transition in the bulk structure, most likely a 
percolation transition at which the surfactant structure changes from a spatial distribution of individual micelles to a 
percolating network. These different structures apparently do not seriously affect the mechanical properties of the gel 
as long as the surfactants are in place. However, once they are soaked out, the resulting percolating void network is 
likely to weaken the gel and lead to large relaxations as the elastic polymer network readjusts its configuration. While 
the percolation threshold is a non-universal quantity and thus depends on the particular system, for site percolation 
on standard three-dimensional lattices it lies roughly in the range from 20 to 30 percent by volume (^?J. This is 
consistent with the range in which the surface structure and the opacity of the templatedgels undergo their strongest 
changes. Previous work on the phase diagrams of alkyltrimethylammonium surfactants Eq] suggests a transition from 
micelle to hexagonal phase at around 35 wt% TTAB at room temperature. 

Further support for the percolation hypothesis is provided by recent computer simulations of a two-dimensional 
frustrated spring network with a free surface and a variable density of holes |p3[ |29f , which were inspired by the 
experimental results presented here. In those simulations the structure of the free surface remains relatively constant, 
with a self-affine structure on short length scales and a finite saturation width on large length scales, as long as the 
volume fraction of holes is far below the percolation limit. As the hole concentration approaches the percolation 
threshold, both the the saturation width and the crossover length increase in a manner qualitatively similar to our 
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experimental observations. The saturation width, in particular, appears to follow the same power-law divergence as 
the bulk correlation length over a wide range of hole concentrations approaching the percolation threshold. These 
observations form our justification for the association of w sa t with the bulk correlation length in the estimate of the 
change in the turbidity, discussed above. 



V. CONCLUSION 



We have obtained AFM images of relaxed, free surfaces of PAAm hydrogels templated with TTAB surfactant and 
performed a numerical scaling analysis of the observed surface structure. In doing so, we have developed methods 
that can be used not only for our data, but also to perform scaling analysis of any surface imaged by an AFM. We 
showed that the calculation of the two-dimensional surface width and one-dimensional correlation functions gave con- 
sistent, useful results. Two-dimensional box-counting and spectral methods were unsatisfactory, while results from 
one-dimensional surface-width measurements provided over-estimates at small length scales. The one-dimensional 
increment correlation function in the fast scan direction, -fff as t( r )i is simple to calculate and gives a reliable char- 
acterization of rough surfaces. By calculating this correlation function we have determined that the surfaces of the 
templated gels are self-affine over a significant range of small and intermediate length scales, with a roughness (Hurst) 
exponent in the range 0.85 to 1. This self-affine regime is terminated by a cross-over length, above which the surface 
width is scale independent. The increases in the cross-over length and the saturation width of the surface with in- 
creasing surfactant concentration coincides with the observation of a large increase in the scattering of visible light. 
We propose that the observed effects are caused by a percolation transition in the bulk surfactant system, from iso- 
lated micelles to a percolating network. Further investigations of the bulk structure of the gels at different surfactant 
concentrations are planned. 
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APPENDIX A 

The theory of light scattering by refractive- index fluctuations leads, in the Born approximation, to the following 
expression for the spectral density of light scattered from a volume V, a distance R away in the far-field zone p0| : 

m=io^^^s^(k) ) (ad 

where 

S {Sn) (k) = { d 3 rG(f)exp[ik-r\ . (A2) 



Here Jo is the intensity of the incident light, A is the wavelength of the scattered light, i\) is the polarization angle, 
k is the scattering wave vector, and G(r) = (Sn(ro)Sn(ro + r)) / ((Sn) 2 ) is the normalized correlation function for the 
refractive- index fluctuations Sn — n — (n) . 

We can calculate the total scattering intensity at a particular wavelength A by using the relation |fc| = sin(|), 
where 9 is the scattering angle, and integrating I(k) over all angles, i.e., 

/(A)=/ ^A- 1 (A,0 (A3) 

where the turbidity or inverse attenuation length is 

hT l (\,£) = J d6d<t>sin6~(l + cos 2 9) S {Sn) (k) , (A4) 
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and £ is the bulk correlation length. The term in the integral of the form (1 + cos 2 0)/2 corresponds to unpolarized 
incident light [ j3l| . Assuming an isotropic exponential correlation function, G(r) = exp [— r/£], we find 



A- x (A,0 



4ir 3 ((6n) 2 ) a 2 



2 2 In (a 2 + 1) 



(A5) 



where a = Aw^/X. Alternatively, using a Gaussian correlation function, G(r) = exp [— (r/£) 2 ], we find 



27T 7 / 2 ((fa) 2 
Ao 1 



{(64- 



+ a 4 ) - (64 + 8a 2 + a 4 ) exp [-a 2 /4]} 



(A6) 



A reasonable value for the mean refractive- index flu ctuati on b etween the PA Am and water regions in a gel is ((6n) 2 
0.01. The resulting turbidity, as given by Eqs. (A5) and (A6), is shown in Fig. Rq. 
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TABLE I: Number of data sets (images) used for analysis. 





Scan Size L [/im] 


% Surfactant 


1.25 


2.5 


5.0 


10.0 


20.0 


Total 


% 


6 


5 


4 


3 


2 


20 


10 % 


5 


5 


5 


3 


5 


23 


20 % 


5 


3 


5 


5 


6 


24 


30 % 


4 


5 


5 


4 


2 


20 


40 % 


4 


5 


4 


4 


2 


19 


Grand Total: 


106 



TABLE II: Summary of results of scaling analyis for templated gel surfaces. 



% Surfactant 


Roughness Exponent a, from: 


RMS Width 
Wsat [nm] 


Cross-over length 
lx [nm] 


Two-dimensional (wj) 


One-dimensional (wi) 




% 


0.41±0.02 


0.77±0.02 


0.63±0.06 


1.3±0.1 


260±70 


10 % 


0.56±0.03 


0.90±0.02 


0.73±0.02 


2.3±0.4 


260±20 


20 % 


0.86±0.02 


1.009±0.006 


0.89±0.02 


7.4±0.3 


260±10 


30 % 


0.80±0.01 


1.023±0.005 


0.87±0.02 


41±3 


400±40 


40 % 


0.86±0.02 


1.026±0.005 


0.83±0.02 


80±10 


540±30 
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FIG. 1: Schematic representation of templated-gel synthesis, (a) In step one, monomers, cross- linkers, and surfactant micelles 
are mixed in solution, (b) In step two, polymerization occurs around the unreactive surfactants, (c) In the final step, the 
surfactants are removed by soaking in distilled water, leaving behind templated voids. The thin dashed lines represent chemical 
bonds between monomers and cross- linkers. The remaining symbols are explained in the figure legend. 




FIG. 2: Schematic representation of the sample holder designed for AFM imaging of hydrogels. 
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FIG. 3: AFM images of a set of successively zoomed-out images centered at the same location on a 20% templated gel surface. 
Parts (a) through (e) represent L x L scan sizes with L = 1.25, 2.5, 5.0, 10.0, and 20.0 (im, respectively. The gray scale 
represents the surface height, with white (black) corresponding to the maximum (minimum) height. The lateral drift of the 
sample was very small. Note the increase in height range with increasing L. The white squares indicate the next smaller scan 
areas. 




FIG. 4: AFM images of a 40% templated gel surface at L = 1.25 /jm scan size, (a) forward scan direction, and (b) reverse 
scan direction. In each image white (black) corresponds to the maximum (minimum) height. The distortion in the direction of 
the tip motion (indicated by white arrows) is readily observed. 
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FIG. 5: AFM images of gel surfaces with surfactant concentrations of (a) 0%, (b) 10%, (c) 20%, (d) 30%, and (e) 40%. The 
scan size is L = 10.0 fim with 512 x 512 pixel resolution. In each image white (black) corresponds to the maximum (minimum) 
height. The height range depends strongly on the surfactant concentration, as is also shown in Figs. |^ and ^| 




FIG. 6: Three-dimensional views of AFM images of gel surfaces with L = 10 /im and surfactant concentrations of (a) 20%, 
and (b) 40% (corresponding to Fig. [| (c) and (e), respectively). The perspective shown is that of an observer sitting at the 
center of the image at a height of 1 /xm above the average height of the surface and looking down at a corner. The light source 
is behind the observer and to the left. 
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FIG. 7: Plots of typical scan lines for L — 10.0 /im gel surfaces (from Fig. g) of various surfactant concentrations. Note the 
strong dependence of the height range on the surfactant concentration. 




FIG. 8: Structure factor S(k x , k y ) of a 40% templated gel on a logarithmic gray scale. The structure factor is given in arbitrary 
units. The anisotropy in the AFM images, which appears as large ridges along the k x — and k y — axes, makes circular 
averaging unreliable. 
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FIG. 9: A typical plot of (wi) vs I on a log- log scale, obtained from two-dimensional images for (a) gel surfaces of various 
surfactant concentrations (corresponding to Fig. H) and (b) 20% templated gel surfaces at various scan sizes L (corresponding 
to Fig. |^). The solid lines in (a) represent linear fits to selected data indicating power-law behavior for limited ranges in I. 
Note the nice overlap of data for various scan sizes in (b). 
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FIG. 10: A typical plot of the one-dimensional flf as t(r) vs r on a log- log scale for (a) gel surfaces of various surfactant 
concentrations at scan size L = 10 /im (corresponding to Fig. ^) and (b) 20% templated gel surfaces at various scan sizes L 
(corresponding to Fig. []). The solid lines in (a) represent linear fits to selected data indicating power-law behavior for limited 
ranges in r. Note the nice overlap of data for various scan sizes in (b). This figure should be compared with Fig. |^. 
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FIG. 11: Plots of the one-dimensional H s \ ov ,(r) and H f ast (r) vs r on a log- log scale for L — 10 fim (a) 0% and (b) 40% 
templated gel surfaces. It is clear that for the untemplated gel surface the signal-to-noise ratio is poor. For the 40% templated 
gel surface, the noise-induced random offsets between scan lines are smaller, but not negligible. 
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FIG. 12: Plot of dlog 10 [f/f ast (r)/2]/d[log 10 r] for L = 10 /im. The inflection point of this curve gives an estimate for the 
lateral cross-over length, i x , which is listed in Table III 
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FIG. 13: Roughness exponent a vs surfactant concentration. The measured values of a are averages over the different scan 
sizes for each surfactant concentration. The error bars represent the standard error. 
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FIG. 14: Saturation surface width w sa t vs surfactant concentration. The error bars represent the standard error. 
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FIG. 15: Cross-over length l x vs surfactant concentration. The error bars represent the standard error. 
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FIG. 16: Turbidity A _1 (A,^) vs correlation length £ for wavelength A = 550 nm in the visible region. The solid and dotted 
curves represent an exponential and a Gaussian form for the correlation function, respectively. Note that the turbidity changes 
by about six orders of magnitude as the correlation length changes from 1 nm (corresponding to w sat for a 0% templated gel) 
to 80 nm (corresponding to w sa t for a 40% templated gel). This increase in turbidity indicates a change from a transparent to 
an opaque material. 



